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Abstract 



We consider large N Yang Mills theory with D adjoint scalar fields in d 
dimensions for d = or 1. We show the existence of a non-trivial saddle point 
of the functional integral at large D which is characterized by a mass gap for 
the adjoint scalars. We integrate out the adjoint scalars in a 1/D expansion 
around the saddle point. In case of one dimension which is regarded as a 
circle, this procedure leads to an effective action for the Wilson line. We find 
an analogue of the confinement /deconfinement transition which consists of 
a second order phase transition from a uniform to a non-uniform eigenvalue 
distribution of the Wilson line, closely followed by a Gross- Witten-Wadia 
transition where a gap develops in the eigenvalue distribution. The phase 
transition can be regarded as a continuation of a Gregory-Laflamme transi- 
tion. Our methods involve large values of the dimensionless 'tHooft coupling. 
The analysis in this paper is quantitatively supported by earlier numerical 
work for D = 9. 
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1 Introduction and Summary 



Matrix models in low dimensions (especially and 1) have served as useful 
tools in many contexts. These include (i) c < 1 matrix models, which cor- 
respond to non-critical string theories [I], (ii) large N reduced models and 
their variants [21 [3], (iii) BFSS matrix theory, which corresponds to DLCQ 
of M-theory [4] , (iv) IKKT matrix model of type IIB string theory [6] , (v) 
DO brane black holes [3 El El EDI [HI H21 [131 E], (vi) KK reduction of 4- 
dimensional M = 4 SYM on S 3 [151 EE El CEE] , (vii) The BMN matrix model 
[T9j . (viii) the matrix model of unstable DO branes [201 12], (ix) D branes 
on tori [22j [231 I2H [25], etc. In many cases, one can regard these models 
as dimensional reduction of large N Yang-Mills theories. In models arising 
from D branes, the YM theories are typically supersymmetric; however, in 
some situations the theory is effectively described by the bosonic sector. The 
obvious advantage of such a description, when it is possible, is that it is easily 
amenable to numerical calculations, and in some fortunate circumstances, to 
some powerful exact methods. 

In this paper, we consider the following system for d = 0,1: 

s^^ld^ilj^D^D^-Y^lr^?], (i.i) 
g j \ i=i i, j j 

where refers to the covariant derivative — i[A^, .]. A^Y 1 are SU(N) 
matrices. For d = there is no gauge field and the first term is absent. For 
d — 1, there is a single gauge field Aq which is non-dynamical. 

The action f ll.lD can be regarded as a dimensional reduction of D + 
d dimensional bosonic YM theory to d dimensions. The specific physical 
context we have in mind in the present paper is related to case (ix) above, 
which is discussed in j22j [23j EU [25] and reviewed below in Section [51 

In the d = 1 case, we consider the dimension as a circle, of circumference 
(3, and study the partition function and other quantities as a function of j3. 
It was conjectured in [22j [231 [2U [25] , on the basis of numerical investigations, 
that the d — 1 system exhibits a phase transition which is analogous to the 
confinement /deconfinement transition of M = 4 SYM on S 3 . The phase 
transition was argued to be the weak coupling continuation of the black 
hole/black string transition of the d = 2 model. One of the motivations 
for the present work was to understand the nature of the phase transition 
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analytically. 

Our results are briefly as follows: 

1. In the limit of large D, the model (11.11) has a non-trivial saddle point 
characterized by a non-zero value of (TyY^ 1 ). The large D scaling is 
defined by keeping a modified 'tHooft coupling A = \D fixed. 

2. In this limit, fluctuations around the above saddle point are suppressed 
by powers of 1/D. This enables us to develop a systematic expansion, 
in 1/D, of the partition function and other quantities as a function of 
the radius of the circle. 

3. In the d = case, the exact partition function is calculable up to 1/D 
at finite N (first shown in [6]), whereas in the d = 1 case our results 
are obtained in the leading large N limit. 

4. Since the condensate provides a dynamically generated mass to the 
adjoint scalars Yj, it is possible to explicitly integrate them out. In the 
d — 1 case, this allows us to compute an effective action S c s(W) for 
the Wilson line 



in a 1/D expansion. 

5. The S e ff computed in this fashion provides the first analytic evidence, 
in a 1/D expansion, for the phase transitions mentioned above. It con- 
firms the appearance of a double transition^: (i) a second order phase 
transition characterized by the onset of non-uniformity in the eigen- 
value distribution of W, followed by (ii) a third order Gross- Witten- 
Wadia (GWW) transition signalling the appearance of a gapped phase 
[261 1271 128] . The appearance of a double transition is supported by 
the numerical works in [2U 125^1 . The phase transition temperatures 
T c \ and T C 2, computed up to 1/D, show excellent agreement with their 
results in the D = 9 CcLSC, clS shown in the following table (see Section 
14.31 for more details and other comparisons): 

1 We thank S. Minwalla for a discussion on this point. 

2 We differ from [24j[25], though, regarding the order of the two transitions. See Section 
E3] for details. 




(1.2) 
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T cl 


T c2 


Our result 


0.895 


0.917 


Numerical result 


0.8761 


0.905 



6. Our methods involve large values of the dimensionless 'tHooft coupling. 

The large D technique used in this paper for d = 0, 1 has earlier been used 
in [6] in the d = context. Our results for d = are in complete agreement 
with those of [6] , though our method is slightly different (see Section [3] for 
details) in a way that enables a natural generalization to higher dimensions. 
A large D expansion has also been used in certain lattice theories [29]. Our 
methods also have some overlap with that of [H [9] where a series of self- 
consistent equations (gap equations) are introduced to determine various 
condensates and a GWW type transition is suggested. 

The paper is organized as follows. In Section |2] we set up the main cal- 
culational method for the d = 1 model. The method consists of introducing 
an S'0(-D)-invariant dynamical field to get rid of the commutator-squared 
interaction. This leads to an action quadratic in the adjoint scalars with a 
dynamical mass term, which allows us to integrate them out. We work out 
the example of the d = model in Section [3] to test the method. We show 
the existence of a non-trivial saddle point in which a mass is generated for 
the adjoints. We explain the nature of the large D limit and compute the 
partition function at finite N in a 1/D expansion. In Section H] we come 
back to the d = 1 model to derive the effective action for the Wilson line in 
a 1/D expansion. We show the existence of a second order phase transition 
followed by a GWW transition. In Section |5] we provide a D brane realization 
of our model, following earlier work [221 I231 I2U |25] which we briefly review. 
In Section [6] we conclude with a discussion. 

This paper arose as part of a larger project of exploring dynamical black 
hole/black string transitions in terms of a dynamical large N transition in a 
unitary matrix model [30] . 

2 The d = 1 model: preliminaries 

We will set up the d = 1 model in this section, test the formalism with the 
simpler case of d = in Section [3j and continue on to Section @] to solve the 
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d = 1 model. A D brane realization of the d = 0, 1 models is discussed in 
Section 

It is convenient to rescale the adjoint scalars Y 1 in (11. ip to gY 1 , so that 
the action becomes 

S=j\ Tr 1 (W) 2 - ^ ^ [Y', y") [Y', Y") \ . (2.1) 

We have assumed the theory to be on a circle, of circumference /3, which can 
be interpreted either as a Euclidean time direction or as a spatial circle. The 
covariant derivative is defined by DqY 1 = d t Y T — i[A , Y ]. 

The large N limit is defined by keeping the 'tHooft coupling A = g 2 N 
fixed, as N — > oo. It is convenient to define the following related dimension- 
less quantities 

/3 eff = /3A 1 / 3 = /3(s 2 A0 1/3 , (2.2) 
A cff = A/3 3 = /3 3 ff . (2.3) 

For later convenience it is useful here to summarise various alternative def- 
initions of coupling constant which we will use in this paper in different 
contexts: 



A 'tHooft coupling g 2 N 

g 2 'tHooft like coupling at large D and finite N g 2 D 

^ 'tHooft like coupling at large D and large iV XD 

A e g dimensionless 'tHooft coupling A/3 3 



The last line refers only to the d — 1 case. 

We would like to explore the partition function 

Z = J VA VY T e- s , (2.4) 

as a function of /3 e g. We will also be interested in the effective action S e s(A ), 
defined by 

exp[-S , eff (A )] = J VY T e- s . (2.5) 

It will turn out that S e g depends only on the gauge-invariant content of the 
gauge potential, namely on the eigenvalues of the Wilson line W (II. 2(1 . 
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The first step in solving the model consists of making the action (12. ip 
quadratic in the Y 1 by introducing an auxiliary field B ab . Let us write 
Y 1 = J2a=i l Ya^a, where A a are the generators of SU(N). This leads to an 
expression 

- Tr\Y' t Y J ]\Y I ,Y J ] = {Y^)M ab}Cd {Y c J Yi), (2.6) 
which is written in terms of S'0(-D)-invariant F-bilinears, where 



M n 



,,cd = -i|Tr[A a , A c ] [At, A d ] + (a <B- b) + (c <B- d) + (a <R- b,c d)|. 



(2.7) 



Properties of the matrix M are discussed in detail in Appendix |A] Using the 
fact that M is invertibldlR, we can write 



Z=M j VBVA^VY 1 e~ s(BMy) , 
S(B, A , Y) =J dt \ (DoY*) 2 - i\B ah YlY* + -L Bab M^ cd B cd 

"(2.8) 

with the following classical equation of motion for 

J 2 Ma b ] cd B cd = iY J a Yl (2.9) 

In the above 1/W = / VB exp[- / dt B ab M^ cd B cd /(4g 2 )}, which we will 
ignore in the rest of the paper since it involves only a numerical factor. Since 
the action (I2.8P is quadratic in Y 1 , we can formally integrate them out, to 



3 It might be puzzling, at first sight, that M does not have a zero mode since 
clearly vanishes for special configurations of the 'Y 1 ^, e.g. when all 'Y 7 's commute. The 
resolution is that M has both positive and negative eigenvalues (see (|A.13p ) and therefore 
has light-like vectors which, however, do not correspond to zero eigenvalues. We thank 
Toby Wiseman for a useful discussion on this point. 

4 The indefinite signature of M (|A.13j) involves some subtlety in choosing the contour 
of the functional integral in (|2.8[) . E.g. we need to choose real contours for components 
of B a i, along the positive eigenspace of M and purely imaginary contours otherwise. The 
correct choice ensures finiteness of the normalization constant Af. 
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get 



VBVA e~ 



-S ea (B,Ao) 



S cS (B, Aq) = [ dt l^B ab M; b \ d B cd + ^ logdet ({D 2 ) ab + iB ab ) . (2.10) 
Jo ^9 * 



For large DUu which scales as 1/g 2 , we may expect the one-loop determinant 
to be comparable with the classical term in (I2.10p and hence modify the naive 
classical solution B ab = 0. We will assume, and shortly justify, that ( I2.10p 
admits a gauge-invariant time-independent solution, of the form (see, e.g. 
Eqns. <KWf and < KE\f ) 

B ab = iA 2 5 ab . (2.11) 

The appearance of i V is due to the fact, as we will see later, that the solution 
corresponds to a saddle point in the complex plane. The condensate in terms 
of the original physical variables Y 1 , however, turns out to be real: 

(TrY'y') = ^A 2 , (2.12) 

where we have used (12. 9p . and also M~ h l cd 5 cd = j^S ab , which is derived in 
Eqn. (tP2|) . 

To proceed, we write the B ab field as the sum of a constant trace piece 
and the rest, as 

B ab (t) = B 5 ab + gb ab (t). (2.13) 

where b ab (t) satisfies J dtb aa (t) = 0. We will show below that B has a saddle 
point value of the form B = «Aq, consistent with ( 12. lip . Substituting ( 12 . 1 3[) 



In much of this paper, we will treat D as an arbitrarily specifiable parameter, except 
in the section dealing with comparison with D branes, where we put D = 9. The precise 
scaling at large D is defined in Eqns. p.6j) and (|3.f 2[) . See also comments at the end of 
Section Ol 

6 Although the large TV scaling is not apparent in (|2.f 0[) , the action as well as the 
measure admits a topological expansion in powers of 1/N for fixed A, as is apparent from 
expressions such as (|3.13p . 
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in the action in ( 12. 8 p we get@ 

Z = JdB VA Vb ab VY 1 e~ s , S = S + S q + S int , (2.14) 

wher^l, 




(2.15) 



Before proceeding to solve this model, let us discuss the d = matrix model 
CIS db partial test of our formalism. 



3 The d = model 

The d = model is defined by a partition functional 

Z = I dY 1 expt-S], S = - 9 \rJ2[y I ,y J ?- (3-1) 

Here Y" 1 are S 1 ?/ (iV) hermitian matrices which are normalized the same way 
as in ( 12.11) . A D brane interpretation of this model is discussed in Sectional 

7 Note the appearance of an effective mass term for the adjoint scalars Yg in the saddle 
point Bo = iAq. In [25 such a mass term is added by hand to integrate out the Y 1 ; here 
the mass term is dynamically generated. 

8 Strictly speaking, the last term in is a cubic term and should be regarded as an 
interaction. However, when we consider Bq as an external parameter (unintcgrated) we 
can regard this term as quadratic. 

9 Most of the results in this section are in [6] who have discussed this model earlier in the 
context of the IKKT matrix model. Our method, however, is slightly different, especially 
in the way we distinguish between the diagonal and the off-diagonal fluctuations of the 
auxiliary field B a f, (see, e.g. (|2.13[) and (|B.2[1 ) which allows for a natural generalization to 
the d = 1 model. 
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Like in the previous section for d = 1, we can rewrite (13.1 j) as 
Z = / dBodbabdY 1 exp[-S], (3.2) 



5* — 5*0 + S q + S'int, 

5o = ^' ^ = l b * M *J>cd ~ ^Yl\ S mt = -jb^Y, 1 . (3.3) 

Integrating out the Y 1 gives us the d = analogue of (I2.10p where B ab is 
split into So and b a b- 

dB db ab e- SMhKb \ 



N B 1 D 

S cf{ (B , b ab ) = -^f + -babM- b ] cd b cd + — log det (iB 5 ab + ^6 a6 ) . (3.4) 



3.1 The large D limit 

Let us make a formal Taylor expansion of S e g in powers of b ab : 
S cS (B , b ab ) = S<f>\B ) + \s% cd (B )b ab b cd + S™\ 



NBn D 



S®{Bo) = —^ + - log det (iB S, 



8g 2 2 



ab 



(2) 1 . # 2 L> ^ log det (iB 6 r8 + ib rs ) 



2 ab ' cd ' 2 960696c 



b rs =0 



S£ = 0(b% (3.5) 

where in defining we have used b a b = gb ab in order to exhibit the g 
dependence explicitly. There is no linear term in the above Taylor expansion 
since b aa = 0. Let us define a large D limit by keeping 

t = g 2 D, (3.6) 

fixed. It is easy to see that 

S^(B ) = O(D), S a f,(B ) = O(l), Sg = 0(l/D). (3.7) 
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Let us now do the integral in (13 .4p over b a b, to give 

Z=J dB exp[-S(B )], (3.8) 

--^(Sb.sO + ^CsD + of-^y (3-9) 



D " K u '" / D \D f 

where 

S l {B (h ~9) = \\ogdetS% d . (3.10) 

The first 'log det' is essentially a 1-loop integral over Y 1 , while the second 
'log det' is a 1-loop integral over b a b (see Figure [6|a) and (b)). In Appendix 
IB1 we will present an explicit computation of these quantities. We find (see 

(E2J) 

NB l , ( N2 ~ l )^J B l 



S (B ) = —± + ± -^log 



8g 2 4 & V 9 2N 



?fR , N^-l ( g 2 N\ N\N + l)(N-3) , ( 2~g 2 \ 

Si{Bo) = —^ log y - -Bi) + § log y 1 - b!) 



+ N2{N -l )iN + 3) lo g (l + %). (3.11) 





In the d — 1 case, an explicit finite N, large D result such as Eqn. (13.111) 
is difficult to obtain, but we will derive an analogue of Eqn. f)3.14p below. 
Furthermore, as remarked at the end of Section 14.11 for d = 1 we will not 
take the strict D = oo limit since criticality involves 1/D effects. 

3.2 Large D, N limit 

It is easy to see that both Sq and Si admit a 'tHooft-like expansion in which 

A = \D = g 2 ND = g 2 N, (3.12) 
is kept fixed. We obtain an expansion of the sort 

5 ,o + -r^So.i + •■■)+ TT (^i,o + -t^Si.i + •••)+•••• ( 3 - 13 ) 



DN 2 V ' N 2 ' D\ ' u N 2 
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In the diagrammatic evaluation described in the Appendices, such an expan- 
sion indeed corresponds to a topological expansion. Explicitly, from (13. lip . 

we get 

S Bl 1 -Bl 1 ( A 1 ( A V 1, / A 

+ O (1/D 2 ) + O (1/N 2 ) . (3.14) 



3.3 The saddle point 

We are left with evaluating (13. 8p . Because of the appearance of an overall 
factor of N 2 in (13.141) . we can evaluate (I3.8P using a saddle point method 
(see a similar calculation presented in Appendix O in a toy example). The 
saddle point value is given by 

B = iAl, A* = 2\(l + 1 L S )+0(1/D 2 ). (3.15) 

The same result is also derived in [6] in a slightly different manner. [6] 
also performed a numerical analysis which agrees with the above analytical 
calculation and also with the numerical calculations of [22| 123] . 
Using the above saddle point, we get the free energy, 

„ ,o g Z 1 log2 + 1, 5 1 3A +0 MV „ 16) 



DN 2 4 4 D \ 8 2 2 J \D 2 , 

4 The d = 1 model: result 

After gaining some experience with the d = matrix model, we now return 
to the more involved case, the d = 1 model. We start with (I2.15p . and as with 
Eqns. (13. 8p and (1B.2j) . we first integrate out the Y 1 and the b a b to obtain 



J dB VA e~ s( - Bo ' Ao \ (4.1) 



where 



e -S(B ,A ) 



J VbabVY'e- 3 , (4.2) 
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with S defined as in (12.151) . Different from the previous section, we consider 
large N case only in this section. 

It is convenient to parametrize Aq by choosing a gauge in which Aq is 
time-independent and is also diagonal: A 0i j = acidy. The gauge-invariant 
content of A Q is then given by the moments 

1 1 N 

Wn = W" = -5>^, (4.3) 

i=l 

where W is the Wilson loop operator, defined in (11. 2p . The above gauge 
fixing gives rise to a Faddeev Popov Jacobian (See [I~6] ) 

VA = 11 da ie ~ SFP , S FP = N 2 J2 l\ u n\ 2 - (4-4) 

i n 

It is convenient to parametrize B = iA 2 , since the saddle point value 
will be real in terms of A as in the d = case (I3.15p . From now on, we will 
denote S(B , A ) as <S(A, {«„}). 

Note that there is a Jacobian involved in changing from the integration 
measure over the eigenvalues «j to the integration measure over u n ,u n ; how- 
ever, it is O(N) and is hence subleading compared to the classical action 
which is 0(N 2 ) [31]. Since in this section we will be concerned with the 
leading term in the 1/N expansion, we will ignore this Jacobian. 

4.1 Computation of <S(A, {u n }) in leading large D 

As in Sectional we can ignore the interaction Si nt in the large D limit. Hence 
the leading result of the effective action is obtained by integrating out the 
Y 1 from S q in f[2~TT5]) . 

We can integrate out Y 1 by using the propagator studied in Appendix [Dl 
(following [12]) and obtain 

D 9 a9XX DN 2 /3A Ai", „ . . 

- log (det {-Dl + A 2 )) = A— -DJ2 -K| 2 - (4.5) 

n=l 

Here x = e _/3A and we have ignored 1/N terms and irrelevant constant terms. 
We also ignored a temperature dependent divergent term. 
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Combining the above equation with So from (I2.15p . and adding the con- 
tribution from ( I4.4p . we get 



<S(A,K}) /3A 4 , 06 ^f l/D-x n \ , |2 

— + ^( K|- (4.6) 



DN 2 8A 2 n 



n=\ 



The 1/-D term comes from (14.411 . The reason it is kept here is that near the 
critical temperature this term will turn out to be more significant than other 
0(1/ D) terms from S,;„ t which we will encounter in the next subsection. 

Our task is to evaluate (14. 1 j) . with S(B ,A ) = <S(A, {u n }) given above. 
It is useful to first perform the integral over Bo, using a saddle point method 
similar to Section [3j In other words, for fixed external u n , let us now solve 
the saddle point equation 



A 3 
2A 



oo 

i + 5>-^ A K| 2 = 0. (4.7) 



n=l 



It is difficult to solve this equation for A({u n }) exactly. However, for small 
u n , it is possible to solve it in a power series in the u n . This leads to the 
following saddle point solution 

A (K» = A 1/3 + f E S >-I 2 j + • • • ' ^ 

where 

x = exp[-/3A 1/3 ]. (4.9) 

Substituting (14. 8 \ in (14.6)1 . we get a Landau-Ginzburg type effective action 
for the u n : 

= I/3A 1 / 3 + ai \ Ul \ 2 + fciH 4 + f» n | 2 + • • • , 



DN 2 

1 



n=2 



(l/D-x n ) 
n 

1 „r 1 



6i = -/3A 1/3 x 2 , (4.10) 
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T = T C 



cl 




S/DN 2 



Figure 1: Phase transitions: S vs (see (I4.10p . ( I4.26jl ). As T crosses T&, 
Ui becomes tachyonic and there is a second order phase transition which 
signals an onset of non- uniformity in the eigenvalue distribution p(a). At 
T = T C 2, characterized by a potential minimum at |iti| = 1/2, a gap develops 
in the eigenvalue distribution, signalling a GWW transition. 

where the • • • involve other u\ terms for n > 1, which are ignored for reasons 
stated below. 

Let us analyze the phase structure of the theory by using (14.101) (see 
Figure ITj)J*°l Our analysis will be similar to [T5]. Note that for x < 1/D all 
'a ra 's are positive. This implies that {u n — V n — 1, 2, ...} is a minimum of 
the potential 

Recall that U\ — = TrVT corresponds to an analog of the confinement 
phase in gauge theory. The vanishing of all u n also has a familiar interpre- 
tation. Let us define an eigenvalue density of the Wilson line (II. 2p by 



10 We will show below that inclusion of higher loop terms does not change the nature 
of phase transitions, although it changes the critical temperature and numerical values of 
various thermodynamical quantities. 

n The issue of whether it is a local or a global minimum is more subtle, and depends on 
details of higher order terms in (|4.10[) . We will argue below that in a 1 /D expansion the 
higher order terms can be ignored and u n = is a global minimum. 




(4.11) 



n=l 
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In terms of this, the u n are given by 

Un = I dap{a)e- inf3a . (4.12) 



The vanishing of all u n therefore corresponds to a uniform eigenvalue distri- 
bution. 

Let us now consider the effect of increasing the temperature, or equiva- 
lently increasing x. As x crosses 1/D, i.e. T crosses a critical temperature 
given by 

T c i = = A 1 / 3 -^-, (4.13) 

log D log -ZJ 

the sign of a\ in f l4.10p flips, while the coefficient of \ui\ A and the coefficients 
of all \u n \ 2 ,n > 1 remains positive. The '|w n |'s remain zero for n > 1 
whereas \u\\ assumes a small non-zero value 



/I \\ f~^~ (3D log D)6T 3 /— n ,5/ 2 /5T\ 5/2 

(4.14) 

at T = T c i + for small and positive ST. E§ 

The order of the phase transition can be determined by studying the free 
energy. It is easy to show that, for small ST, the Landau Ginzburg free 
energy is of the form 

S/(DN 2 ) = constant + constant (5T) 2 6(5T) + O -^j . (4.15) 

The second derivative of the above function with respect to ST is discon- 
tinuous at ST = 0. Thus we have a second order phase transition. We 
should remark that the transition is characteristic of the large iV limit and 
is expected to be smoothened at finite N, as has been argued in [16] . 



12 Bccause of couplings such as U-2u\ discussed in Section [4.21 higher it n 's pick up some 
non-zero values at higher orders in 1/D. However, these can be ignored in the present 
discussion. 

13 ST here is assumed to be small enough such that (\ui\) satisfies the bound < 1/2 
discussed below. 
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As we increase the temperature further, we encounter another phase tran- 
sition. To understand this phase transition, we first note that when u n = 
for all n > 1 (which we expect to hold as l ong as x does not cross 1/yD), 
the eigenvalue density can be represented as 14 l 

p(a) — — (1 + 2|mx| cos(/3a)) . 
2n 

As \ui\ increases from small values to 1/2, the eigenvalue density vanishes at 
(3a = Ti. In the present case, the saddle point value (I4.14p . reaches 

the value 1/2 (see Fig. [TJ) when T equals a critical temperature 

A 1 / 3 

T c2 = T cl + . (4.16) 

6D log D y 1 

As T crosses T c2 we have a GWW type phase transition [2^ [27] from a gapless 
eigenvalue distribution to a gapped one. The nature of this transition has 
been discussed in detail in [161 IUJ> where a Landau- Ginzburg potential of 
the form ( I4.10p was assumed, with vanishing u n ,n > 1. Our analysis above 
supports this assumption, hence we can use their analysis. Following Eqn. 
(6.18) in [IB], we find that for temperatures just above T c2 , |iti| grows as 



By comparing this equation with the form of below T C 2 (I4.14p . we find 



that the second derivative of |iii| 2 (or equivalently the third derivative of the 
Landau Ginzburg free energy) with respect to temperature is discontinuous 
at T C 2, although the first derivative or the value of |wi| 2 is continuous. To be 
precise, we fincEl 

I |2 1 x 3 m n T ~ T c2 9D(logD) 3 ( T-T c2 \ 2 

lu-i \ — — — D log D — — ~ + • • ■ , 1 < 1 C 2 

I 11 4 2 & A x /3 4 V A 1 / 3 / 

,9 1 3_ ^T-T c2 27D 4 \ogD [T-T c2 \ 2 m m 

kii = - + -DlogD — = — —I— — - +•••, T>T c2 

II 4 2 & A x /3 2 I \ 1/3 J 



(4.18) 



14 By a suitable shift of the origin of the angle a to absorb the phase of u\. 



15 The rate of change of \u\ | 2 below T C 2 is given by the expansion parameter (log D) 2 \T — 
T c2 \/~\ 1/3 while above T c2 it is given by the expansion parameter D 3 \T — T^l/A 1 / 3 . Hence 
|ui| 2 changes at a much faster rate above T C 2- 
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T 

Figure 2: Phase transitions: \ui\ vs T. As T crosses T cl , \ui\ starts growing 
from zero and equals 1/2 at T = T c2 . The transition at T c2 is a third order 
GWW type transition. Because of the sharp change in the second derivative 
of across T c2 within a very small range of temperature (see Eqn. ( 14 . 1 8 j) 
and also footnote [15]) . it almost appears like a discontinuity in the slope of 
the plot. However, when we zoom into a small temperature interval around 
T c2 , the d\ui\/dT is seen to be continuous, as is analytically evident from 
f )4.18p . Although we have plotted vs T here, to facilitate comparison 
with [21], a plot of |wi| 2 vs T shows exactly similar features. 



where we have ignored corrections of 0(1/D 2 ) and terms proportional to 
(T — T c2 ) 3 and higher. Thus, the phase transition at T = T c2 is third order, 
as in the original GWW transition. (See Figure. [21) 

Beyond T c2 it is in principle possible to have further phase transitions to 
multiple-gap phases. In case of unitary matrix models with general single 
trace actions of the form J2 n C n u n+ C -C- it wa s shown in [321 [33] that additional 
gaps can open up in the eigenvalue distribution as the temperature is varied. 
We will not attempt to study this issue in this paper, since the analysis in 
the relevant ranges of temperatures is complicated. The numerical analysis 
of [24], appears to suggest, however, that the only phase transitions in the 
model are the two already discussed above. 



High temperature: Once the temperature becomes very high, /3 c g 1, 
the model again admits analytic treatment. In this region, the eigenvalue 
density approaches a delta function and we can approximate u n — 1. The 
saddle point equation ( 14. 71) then becomes 

A 3 1 e-? A 

-=^ + - + — ^r = 0. (4.19) 

2A 2 l-e~? A V ; 
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We can approximately solve it for small (3 as, 

A=^) . (4,0) 

This result is consistent with the d — condensate in (13 . 1 5f) by identifying 

A| rf=0 as A//3. 

We end this subsection with a few comments: 

(a) Small \u±\ approximation: Since the second phase transition hap- 
pens when \ui\ reaches 1/2, one may worry about the validity of small u n 
approximation in (14.81) . There is no problem, however, since terms involv- 
ing u n in the saddle point equation (14.71) are also suppressed by x n ~ 1/D n 
around the critical temperature. Therefore, even if is not small, our 
analysis is valid. 

(b) Large D limit \s 1/D expansion: It is clear from the phase transi- 
tion temperatures (14.131) and (14.161) that the phase transitions disappear in 
the large D limit. Hence we should simply regard D as large but not take 
the strict D — > oo limit if we want to explore criticality. 

4.2 1/D correction to the effective action 

In the previous section, we have considered the effective theory (14.61) including 
the 1/D term (14.41) from the gauge fixing and discussed the phase transition. 
However, in addition to this 1/D term, other 1/D corrections can arise from 
Smt in (12.151) . It corresponds to S^o in (I3.13P in the dimensional model. 
Hence we have to evaluate them and show that these corrections are sub- 
dominant around the critical point. 

The terms we should look for at this order, in so far as the issue of phase 
transitions is concerned, £11*6 cLS follows. Besides the explicit corrections to 
|tti| 2 and |ui| 4 in (14.101) . the corrections to the gauge-field independent terms 
are also relevant, since they contribute to the saddle point equation (14. 7p . 
Interaction terms like u\u_2 also affect the Landau- Ginzburg type potential 
(14.10p by generating an effective |wi| 4 term. However as we show below 
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Eqn. (lE.12p . the corrections to the coefficient of |wi| 4 from these interactions 
are order 1/D 4 and we ignore them hera^l. Thus the relevant terms in the 
effective action are 



S(A, {u n })/(DN 2 ) = Co + C 2 \ Ul \ 2 + C 4 K| 4 

and we can explicitly calculate them, 

RA 4 RA 
8A -i 
0A 17. 



+ 



D 



2 

1 + 



4A 3 



A 

4A 3 



1 
I 



A 
4"A 3 



(4.21) 



(4.22) 



C, 




(4.23) 



PA 
AD 



4A 3 



A 

4A 3 



4A 3 



x 



A 

4A 3 



A 



4A 3 



x 2 + 0(x 3 ). (4.24) 



Here we have omitted higher x terms, since we are interested in a range of 
temperatures below or around the critical temperature x ~ 1/D details of 
the derivation are shown in Appendix lEl 



16 Although the cubic interaction u\u-2 merely renormalizes the coefficient of the |ui| 4 
term in the Landau-Ginzburg potential (|4. 10[) . one needs to be careful about integrating 
out the U2 consistent with the positivity constraint of p{a). 

17 Note that the small x expansion is valid for low temperatures. It means that our 
analysis works well for large effective coupling A G ff (|2.3j) (as long as it does not scale with 
D). This assumption is in particular valid around the phase transitions. 
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As in f )4.8p . we solve the saddle point equation for A in powers of U\, to 
obtain: 

A 1 ( 7V5 9 \ 2 „ , 

1 + - ^7T-— +-x| Ml | 2 + ---. (4.25) 



A 1 / 3 D \ 30 32 J 3 

Here the |wi| 4 and higher order terms do not affect (I4.26[) and are dropped. 
Substituting this in (I4.2ip . we get 



with 



3 




81 


° 8 


•H 


~64 4 


, 1 

fll =D 


— x — 


A 1/3 /3 


D 



S/(DN 2 ) = f\ 1/3 e + a'^ml 2 + b[\ Ul \ 4 + • • • , (4.26) 

(4.27) 

x, (4.28) 



y J^l- X ._~^l(~ x „ + 391V5 _ 3181 V 2 (4 29) 

1 3 D \ y I 9 300 y 1800 2400 / v ' ; 

It is obvious that these equations constitute 0(1/ D) fractional corrections 
to various quantities appearing in (I4.10p . 

As argued in the previous subsection, the phase transition temperature 
T cl is characterized by vanishing of a' x and T c2 is given by 6' x = —2a[. This 
gives us the following corrected values of the transition temperatures: 

^^ logD ( 1+ i(m-4)\. ( ,30) 



= log 




4( 


203 v 7 ^ 
v 160 ~ T~ 


- /3 c iA 1/3 






logD 


1 




/ 499073 


D 


-6 + 




I 460800 



203^ , „ 1127^ 85051 

log D 1 

480 j B 1800 76800 

(4.31) 

Although the analysis in this subsection leads to subleading corrections to 
the phase transition temperatures, it is easy to see that the nature of the 
phase transitions derived in the previous subsection remains unaltered. 
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4.3 1/D expansion vs. numerical calculation 

In this section, we evaluate the critical temperatures and a few other quanti- 
ties, using results in the previous subsections in the large D expansion, and 
compare them with the numerical results for D — 9, which were studied in 



The d = 1 model was numerically analyzed in [221 123] ( see Section [5] for 
connection to D branes) where it was suggested that the system undergoes a 
weakly first order Gross- Witten-Wadia type transition (characterized by the 
development of a gap in the eigenvalue distribution p(a)). A more detailed 
numerical study [23] subsequently claimed that in stead of a single first order 
transition, it consists of two higher order phase transitions: (a) from uniform 
to non-uniform p(a), followed closely by (b) a GWW type transition in which 
a gap appears. This is in agreement with the picture of the two transitions we 
derived in the previous subsections (see Figure [3]). Let us compare between 
our results and those of [23] in some detail. 

We first compare the two critical temperatures derived from the numerical 
analysis in [23] with our large D expansion. In order to do it, we note that 
the dimensionless temperature defined in [221 I231 [231 |25] and (I2.2p is 

1 1 D 1 / 3 

Tcff = U = AV^ = AW (432) 

In the units A = g 2 N = 1, used by [23], T e s is simply T. By employing the 
same unit, we obtain the critical temperatures as in Table [TJ The leading 
order results in the 1/D expansion are from Eqns. (14. 13ft and ( I4.16jl . and the 
next order is from Eqns. (I4.30p and (I4.3ip . 

Similarly, we can also compare the value of the condensate A 2 and the 
free energy in the confinement phase, which are given by 

R 2 = ^(TtY i Y i )\^ oo = ±A 2 \^ O0 , F = --g^3logZ|^ oo . (4.33) 

Those can be derived from ( 14.81) and (I4.10p in the leading order, and (14 . 2 5 j) 
and (I4.26P in the next order. The results are also summarized in table [U and 
our results agree with the numerical result remarkably welf^l. 



18 Note that we call the critical temperature from uniform to non-uniform distribution 
as T c i and the next GWW type as T c2 . However, in [24], they used the opposite notation. 
R 2 and Fq in (|4.33[) are defined as r 2 , and eg in [24] - 
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T cl 


T c2 


K 2 


F 


Numerical result 


0.8761 


0.905 


2.291 


6.695 


Leading large D result 


0.947 


0.964 


2.16 


7.02 


Large D including 1/D effect 


0.895 


0.917 


2.28 


6.72 



Table 1: Comparison with numerical results derived in [21] for D — 9 and 
our large D analysis. Here we have used A = g 2 N = 1 units. We list the 
critical temperature T c \ and T C 2, and the condensation and the free energy at 
the confinement phase defined in (I4.33p . The first line is the numerical result. 
The values in the second line are the leading large D results. The third line is 
the result including the first 1/D correction. The fractional differences from 
the numerical results can be checked to be order 1/D in the second line and 
1/D 2 in the third line, as expected. 



We note here that although we find excellent quantitative agreement with 
|24j . the more qualitative inferences in [21] regarding the order of the phase 
transitions are different from ours. The phase transitions at T cl and T c2 are 
claimed in [21] to be of 3rd order and 2nd order, respectively, while in our 
analysis they are the other way around. We believe that this difference may 
be due to the fact that in numerical work it is not easy to ascertain the order 
of a transition except when it is a strong first order transition. The transition 
at T cl in our analysis is described by a classic Landau Ginzburg potential 
which describes a second order transition. For a LG potential involving only 
u\ to describe a third order transition as suggested in [21], we need a u\ term 
which is disallowed by the symmetries of the theory. Likewise, a second order 
transition at T c2 is inferred in [21] by noting a jump in d\ui\/dT. We find, on 
the other hand, that there is a sharp, but continuous change in this quantity 
(see Figure Eland Eqn. (14. 18f) ) . We note that our analysis, of course, ignores 
corrections of order 1/D 2 ; however, we do not expect any qualitative changes 
in the above conclusions for large values of D such as D = 9/^1 . 

19 J. Nishimura informed us that he agrees with the conclusion obtained in this paper 
and that the numerical data in [24) around T c \ arc also consistent with a second order 
phase transition. He also mentioned that fitting the data with that assumption leads to 
a slight increase in their estimate on T c i, which further improves the agreement with our 
value of T cl . The estimate of T c2 in [21], on the other hand, does not depend on the 
assumed order of transition in their analysis. Thus, the agreement with our value of T c2 
is not affected. We thank J. Nishimura for providing us with the results of this reanalysis. 
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1/D expansion for small D: Ref. [25] also numerically analysed the 
transition for D = 2 and 3, and found two transitions as in the D = 9 case 
in [23]. In our study too, we have two transitions for D = 2 and 3, since b[ 
in (I4.29P is positive even for these values of D. We summarise these results 
as follows: 





T cl (D = 2) 


T c2 (D = 2) 


T cl (D = 3) 


T c2 (D = 3) 


Our result 


1.4 


1.6 


1.1 


1.2 


Numerical result 


1.12 


1.3 


0.93 


1.1 



The critical temperatures in their numerical results are close to our results 
Our phase transition temperatures also show agreement with the critical 
temperature numerically evaluated in Ref. [23] for D = 4. 

These results suggest that our analysis seems to work even for such small 
values of D. However, we do not have a detailed understanding of such an 
unexpected agreement. 



5 D brane interpretation 

The d = 1 model ( 12.11) appears in various contexts, as mentioned in the 
Introduction. The context closest to the contents of this paper is that of 
[22l |23| 124"] . which we briefly review in this section. 

Let us consider thermal DO branes in R 8 x S 1 . The distribution of the 
branes is dynamically determined and for a certain parameter region, the 
geometry becomes a black string winding around the S 1 . If we increase the 
radius of the S 1 beyond a critical radius, the Gregory-Laflamme instability 
mode [3l] appears and a black hole solution, localized on the S l , is favoured. 
It is argued in [22] (see also [35]) that this black string/black hole transition 
is first order. 

Through gauge/gravity duality [7], we expect this transition to be repro- 
duced by a d = 1 thermal SYM with a compact adjoint scalar at strong 
'tHooft coupling. By using a T-duality [SB], this model can be mapped to a 

20 Their interpretation of the order of the phase transitions is the same as in [21], and 
is different from ours. The explanation of this discrepancy is similar to the D = 9 case as 
mentioned above. 
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■uniform 



Od YM """ \ gapped ^ 
t' \ A' = 1/t 



Figure 3: Phase diagram of the d = 2 SYM ([53]) . Below A' = t', the 
temporal KK modes can be ignored and below A' = 1/t', the spatial KK 
modes can be ignored. Thus the effective d = 1 (bosonic) YM description 
is valid on the right of the curve A' = t' 3 . The overlap of this region with 
the region below A' = 1/t' additionally admits an effective d = description. 
The two phase transition lines below A' = t' 3 are given by X't' = and 
X't' = f/T c 3 2 , where T c i,c2 are given in (14.301) and (14. 3 f h . A similar phase 
structure was earlier inferred in [21] on the basis of numerical analysis. 



2d SYM on T 2 

S=\ fdx fit Tr [-¥%, + \ V D^YjD^Yj - - V[Y 7 , YjA + fermions. 

92 Jo Jo \4 Z i=i 4 i j J 

(5.1) 

This theory is characterized by two independent dimensionless constants: (a) 
A' = \2L 2 where A2 = g%N is the 'tHooft coupling, and (b) t' = L/(3 2 , the 
dimensionless temperature. However, the analysis of this theory at strong 
'tHooft coupling is difficult. Instead of investigating the above transition at 
strong coupling, the gauge theory allows us to study a continuation of the 
phase transition to weak 'tHooft coupling. 

It has been argued in [2"2"| |2"3] that in the range of temperatures given 
by A' 1 / 3 < t' < 1/X' all fermionic modes as well as both the spatial and 
temporal KK modes can be ignored. The theory is then governed by just the 
zero modes which describe the d = model studied in Section |3l (See Figure 
E]). As X't' grows to order unity, the spatial KK modes cannot be ignored 
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any more, though the temporal KK modes can still be ignored. In fact, in 
the range of temperatures and coupling t' 3 > A', the theory (15.1 j) reduces to 
the d = 1 model (12.11) with the spatial circle of length L identified with j3 
of (12 . 1 p . The d = 1 'tHooft coupling A is identified with A 2 //3' so that A c g 
appearing in (12. 3 p is identified as 

A cff = \ 2 L 3 /(3 2 = X't'. (5.2) 

Note that the transitions in the d = 1 model, which we studied in Section 
lU happen around A c fr = w 1,4 from Table [IJ These transitions can 

indeed be regarded as transitions in the d = 2 model (15.11) if t' 4 > A e ff. Thus 
we can reliably expect these transitions to be the continuation of the black 
string/black hole transition to weak 'tHooft coupling^]. 

As has been suggested first in [37], the eigenvalue distribution of the 
Wilson loop (14. lip is related to the geometry of the DO branes. A uniform 
(non-uniform) gapless eigenvalue distribution corresponds to a uniform (non- 
uniform) black string winding around the S 1 , whereas a gapped distribution 
corresponds to a black hole localized on the S 1 . Now, we found in Section [4] 
that the uniform eigenvalue distribution is favoured at low temperature and 
a gapped distribution is favoured at high temperature while a non-uniform 
distribution exists between those two phases. Since the temperature in the 
d — 1 model is mapped to the radius of the original S 1 , there should be a 
phase transition from a black string to a black hole as the radius of the circle 
is increased, which indeed is the case. The fact that our transition consists of 
two closely spaced transitions disagrees, however, with the single first order 
transition in the gravity description. A plausible resolution is as follows. It is 
easy to see that if b\ in (I4.10p is negative when a\ vanishes, there is only one, 
first order, phase transition instead of the two transitions [16]. Therefore the 
gravity picture can be reconciled with the gauge theory calculation if the sign 
of b\ in (14.1 Op flips at some higher value of coupling in the two dimensional 
model. At such a value the two phase transition lines found at weak coupling 
will merge and yield a single first order transition line (see Figure |3]). 

21 The "weak" coupling here refers to the d = 2 coupling A' which satisfies A' < i' 3 (see 
Figure |3]) . We should remark that the analysis in this paper is valid even for large values 
of the d = 1 'tHooft coupling A e fi as we explained in footnote [17] The equivalence with 
the d = 2 model, however, is valid only for temperatures t' A > A c ff . 
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6 Conclusion 



In this paper we have developed a technique of solving matrix models (d = 
0, 1) which are dimensional reductions of D + d dimensional bosonic YM 
theory to d dimensions. The technique involves working in a 1/D expan- 
sion, which allows us to analytically compute free energies and other ther- 
modynamic quantities. In the d — 1 case our results show that the system 
undergoes a double phase transition: a second order phase transition which 
signals onset of a non-uniformity of the eigenvalue distribution p(a) of the 
Wilson line, followed by a third order GWW phase transition which signals 
development of a gap in p(a). Following the arguments in [221 1231 12H |25], we 
interpreted this double transition as a continuation of the Gregory-Laflamme 
(black string/black hole) phase transition to weak coupling. Our results agree 
with the numerical results of [23], and offers an analytic resolution of the issue 
of the order of the phase transitions. 

The large D technique developed in this paper is in principle applicable 
to a variety of bosonic matrix models involving commutator-squared inter- 
actions. The applicability of our techniques would be greatly enhanced if we 
are able to extend them to higher dimensional models d > 1 and to include 
supersymmetry. The extension to higher dimensional models appears to have 
the technical hurdle of computing logdet(D^ + iB ) with dynamical gauge 
fields without making a coupling constant expansion. One possibility is to 
find regions of parameter space in which an effective d — 1 description arises 
(as in the d = 2 toroidal model described in Section |5]) and work around that 
limit. The supersymmetric extension of the large D methods appears more 
challenging, even qualitatively, since the number of bosons and fermions grow 
at different rates as D grows large. We hope to come back to some of these 
issues in a future publication. 

In this paper we have been concerned with thermodynamic properties of 
the matrix models. Another possible application of our methods could be to 
address dynamical questions. Indeed, one of the motivations of this paper 
was to apply these techniques to derive an effective action for gauge fields in 
the time-dependent context and to study dynamical phase transitions using 
this effective action. Work in this direction is in progress [30] . 
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A Some results involving M a ^ c d 

In this section, we will calculate several quantities involving M ab ^ for exam- 
ple, TrM n , M _1 , the eigenvalues of M etc. which are important in solving 
both the d = and the d = 1 models. We begin by first investigating 
algebraic properties of M abtCd . 

A.l Algebraic properties of M a ^ C( i 
As in (12. 7p . M a ^ c d is defined as 

M abM = --{Tr[A , A c ][A 6 , A d ] + (a O b) + (c <-> d) + (a <-> 6, c <-> 

(A.l) 

Mah^d has four adjoint indices and is symmetric under the interchanges a <H> b 
and c -H- d. Hence, in the SU(N) case, we can regard M ah cd as an N 2 {N 2 — 
l)/2 x N 2 (N 2 — l)/2 matrix by identifying ab and cd as two single indices. 
Equivalently, we can regard M as an operator acting on the N 2 (N 2 — l)/2 
dimensional vector space Vb (whose elements can be regarded as B ab ) labeled 
by a symmetric pair of adjoint indices ab, on to the same vector space {i.e. 
M is an endomorphism of the vector space Vb). 
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N-1 



N-1 



= 1 + 



N-1 



+ 



N-2 



+ 



N-1 



N-1 



Figure 4: Irreducible decomposition of the N 2 (N 2 — l)/2 dimensional vector 
space Vb labelled by a symmetric pair of adjoint indices ab, regarded as a 
symmetric product of 'adjoint' x 'adjoint'. The dimensions of the representa- 
tions in the RHS are 1 , iV 2 - 1 , iV 2 (N + 1 ) (N - 3) /4 and N 2 (N - 1 ) (N + 3) /4 
respectively. These correspond to the space obtained from the projection 
through K , Ki, K 2 and K 3 . Note that SU{2) and 577(3) are exceptional. 
In SU(2) the second and the third representations are absent and in SU(3) 
the third one is absent |38|. 



To proceed, let us decompose this N 2 (N 2 — 1)/2 dimensional vector space 
Vb into irreducible representations of SU(N). The decomposition is shown 
in Figure H] and we obtain four irreducible representations. Let us define a 
vector in this space as B a b, then the first 1 in Figure |4]is the 'trace' part B aa 
I 22 ! and the latter three constitute irreducible decomposition of the symmetric 
'traceless' vector. Correspondingly we can define four projection operators 
Kiab,cd {i = 0,1,2,3) acting on this vector space such that K iab cd B cd belongs 
to the i-th irreducible representation, 

B a b = (K B) ab + {KiBU + (K 2 B) ab + (K 3 B) ab . 

We will show that the endomorphism M abfi d can also be decomposed following 
the above equation and hence can be represented as a linear combination of 
the Ki (see ([PPjl ). 

It is possible to construct the Ki explicitly. Let us define the following 

22 Here 'trace' of a vector B ab means the sum of the two adjoint indices, e.g. B aa . 
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four matrices [B]: 

F ab ,cd =^ [Tr (\ a \ b \ c \ d ) + (a^b) + (c^d) + (a^b, cod)], 
G afe , cd =- [Tr (A a A c A b A d ) + (a o 6)] , 

H a b,cd =VabVcd-, 

Iab,cd =2 (<W&d + ^aAc) • (A. 2) 

Those matrices satisfy the following relations, 

2 V A^V 2 V iV 2 / N N 2 

FH = Ml - —z)H, G 2 = I + -^rH - —F, GH = — —H, 

v M 1 ' iV 2 iV AT ' 

H 2 = (N 2 - l)H, FI = F, GI = G, HI = H ,I 2 = 1 (A.3) 

where the product AB is defined by (AB) ahe j = A abjCd B cde f. Because of the 
cyclic property of the trace, this product satisfies AB = BA. Note that I a b,cd 
plays the role of the identity matrix. In order to derive these relations, we 
have employed the identity 

N 2 -l ^ 
a=l 

where i,j are fundamental indices. 

By using these matrices, the four projection operators are represented as: 

2N r F-ln 



K Q = 


N 2 - 






K 2 = 




(' 


- G - 


2 ( 


i' 


N-2\ 


K 3 = 




{' 


+ G- 


2 ( 


i' 


N + 2 \ 



N 2 -A\ N 
F-—L 



N J N(N-l)' 

F-—H) -H). (A.5) 

N J N(N + l) J V ; 

We can show that these matrices satisfy, 

K i K j = SijKi, (A.6) 
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and 



' ab,cd 



Koab,cd + Ki abfid + K 2a b,cd + K 3abcd . (A. 7) 



Thus the K t are indeed projection operators. Actually, we can find the cor- 
respondence between Ki defined in flA.5j) and the irreducible representations 
in Figure HI For example, Ko ab)Cd acts on a vector B cd as 

Koab,cdB c d = jy 2 _ l fiabBcc- (A-8) 

Thus Kq maps the vector to the 'trace' B cc . Therefore Kq is the projection 
operator corresponding to 1 in the RHS of Figure HI Similarly we can find 
the correspondence for other K^. 

In addition to such explicit identifications, we can calculate the traces of 
Ki aS 



Koab,ab — 1, Klab,ab ~ N — 1, 

N*(N + l)(N-3) N 2 (N-l)(N + 3) 

J^2ab,ab = T , J^3ab,ab = ^ • l A - y J 

The values of the traces are equivalent to the dimensions of the irreducible 
representations shown in Figure HI This is another evidence for the corre- 
spondence between Ki and those representations. 

By using (lA.lj) . (1A.2|) and (1A.5[) . M abfid can be described as 

M a b t cd = 2{F ab)Cd — G a b,cd) 

= 2NK 0ab , cd + NK labjCd + 2K 2ab>cd - 2K 3abM . (A. 10) 



A. 2 Results involving 

In this subsection we use the tools developed above to calculate several quan- 
tities associated with M abtCd necessary in the study of the matrix models in 
this paper. 

23 Here 'trace' of an endomorphism matrix A a \y tCd is defined as TrA = A^^i, 
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The inverse M 1 : Equation (12.81) assumes the existence of M 1 . From 
flA.7j) and (lA.lOp . we explicitly obtain 

j 1 1 11 

^ab,cd = l^~jyKoab,cd + J^K labcd + -K 2abcd ~ -^K 3abcd . (A- 11) 

By using this, we can calculate M~ b \ c , which is necessary to derive (I2.12p 
and (12.151) . First we can show Ki abtCC = for i ^ 0, since the irreducible 
representations in Figure H] are 'traceless' except 1. Then we obtain 

M ab,cc = i^K 0ab . cc = 7^5 ab . (A. 12) 

Eigenvalue of M abjCd : We now derive the eigenvector and eigenvalue of 
M a b,cd- By using (lA.lOj) . we obtain the eigenvector as 



M abM {(K B) cd + {KtBU + (K 2 B) cd + (K 3 B) cd } 
=2N(K B) ab + N(KiB) a t, + 2{K 2 B) ab - 2(K 3 B) ab . (A.13) 

where B ab is a general vector. Note that the eigenvalue of (K 3 B) ab is nega- 
tive, which makes the action (12. 8p not positive definite. As we have remarked 
in footnote HI this can be dealt with by making appropriate choices for inte- 
gration contours for different elements of B ab . 

Calculation of TrM n : We now calculate TrM n = M™ 6 ab , which appears in 
the loop calculation of the matrix model. By using Eqns. ( 1A.6I) and (lA.lOp . 
we obtain 

TrM n = (2A0 n TrA' o + JSTTrKx + 2 n TrA^ + (-2) n TrA' 3 . (A.14) 
In the large A" limit, we obtain by using ( 1A.9j) 



TrM = -N-\ TrM = 3N , TrM n = N n+Z (n > 3). (A.15) 

It is also possible to calculate the effective action for finite A" in the d = 
model. To do this, note that in the loop diagrams including b ab , the vector b ab 
satisfies the 'traceless' condition b aa = as in (I2.13p . This condition changes 
the propagator for b ab with M replaced by M'\ 

M' = NK labtCd + 2K 2abM - 2K 3ab<cd . (A. 16) 
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(b ab b cd ) (YlYi) -iba^M 

Figure 5: Feynman rules for the matrix model ( 13.31) . 

Hence the 6 a b-loops actually involve TrM n which are given by 
Mf ■ =N"(N* - 1) + y^ + ^-a) + { _ T ^(N-l) { N + 3 ) 

(A.17) 

Note that it gives the same value to (1A.15I) under the large N limit and the 
traceless condition indeed affects only finite N correction. 

B Details of the d = model 

In this Appendix, we will present details pertaining to the d = model of 
Section [31 Specifically, we will derive (13. lip . 

We start with (13. 3p . The Feynman rules of (13.31) are shown in Figure [5j 
The propagators are given by 

(YX J ) = -^5 ab 5 IJ , (b ab b cd ) = 2M' ab , cd . (B.l) 

The matrix M' ab cd is defined in (1A. 16j) . which is the propagator for b ab satisfy- 
ing the traceless condition b aa = 0. Note that M' abcd is obtained by removing 
K from Mah^d, where K is the projection operator corresponding to the 
trace part b aa (See eq.( lA.8P ). As remarked below fl A. 1 7|) . the difference be- 
tween M and M' appears only at subleading orders in 1/N. 
The effective action S(Bq) in (I3.8P is formally given by 



exp[-S(B )} = J dY 1 dbe- s °- s " U + ^T 



Sint) n 



n=l 



e -S 



iBo) ^ l + DM"^ ■ (R2) 



n=l 



32 



(a) - D 



9 

(b) ~g 2 D~g 2 



9 

(c) ~g 4 D^gyD 



Figure 6: Some examples of vacuum diagrams of (IB. 21) . Figure (a) is the 
leading order in the D expansion (0(D)). (b) is next order (0(1)) and (c) 
is 0(1/D). The diagrams in Figure M also contribute to 0(1). 

where we have dropped some irrelevant normalization factor. Thus, S(B ) 
is given by a sum of all connected vacuum diagrams represented by the 
above equation, with external Bq. The term in the above equation obtained 
by putting S int = corresponds to S^^Bq) in (13. 5p . Diagrammatically it 
corresponds to the y 7 -loop in Figure [6^a) plus terms independent of Y and 
b. Since each Y-\oop gives rise to a factor D, this term is of order 0(D), 
consistent with the arguments in Section 13.11 
Therefore, we obtain S in ( 13 . 1 j) as, 



where g 2 = g 2 D, as defined in (13. 6p . Here we have used the notation 



in anticipation of the fact that the saddle point value (I3.15P will be given in 
terms of real A. 

B.l Calculation of S\: the 1/D correction 

Now we calculate Si in ( I3.10p . 

We can obtain Si directly by using the full propagator for (bb), as in 
Section 13.11 However we calculate it diagrammatically here^l. since this 

24 The connection with Section 13.11 can be made by the formal Schwinger-Dyson sum 

5 (2) = M + AIGAI + MGMGM H where S* (2) is the quantity appearing in (pHij) . while 

G = G( 2 ) appears in (|B.5p . 




(B.3) 



B = iA 2 , 



(B.4) 
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Figure 7: One-loop correction to the (bb) propagators. The correction to (bb) 
is order g 2 . The 1PI propagator comes from this diagram only at this order 
in the 1/D expansion. 




two loop three loop higher loop 



Figure 8: The 0(1) corrections to the effective action in the large D expan- 
sion. 



derivation will be more convenient. It is easy to show that at leading order 
in 1/D, the relevant correction to the 1PI (bb) propagator comes entirely 
from the one-loop diagram of Figure [7J Hence only diagrams described in 
Figure [8] contribute to the effective action. Note that each diagram includes 
planar and non-planar structures in the N counting. Higher loop terms are 
characterized by higher powers of the dimensionless quantity A/ A 4 . However, 
as evident from (I3.15p . this quantity is order 1. Hence we must sum over all 
loops, which we describe below. 

In each diagram in Figure El the two Y ! s in a vertex b^Y^Y^ are con- 
tracted with two Y J s in a different vertex b c dY c J Y/ . Therefore it is convenient 
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to define a composite propagator 

'J\ + (y i y^\Iy, i y j \\ = 



on 

G (2)ab , cd = ((Y*Y c J )(YtYl) + (YttKYiYf)) = --I ab ^ (B.5) 



y..y 

where / a fe,cd is defined in (|A.2[) and it satisfies M' abcd I cd ^ e f = M' abe ^. This 
composite propagator corresponds to one double line loop in Figure and |HJ 
By using this propagator, we can calculate the (n + l)-th loop correction to 
N 2 S l as 

2rj 



(2n 



1 / ■ \ zn 

(02n\ — _ ( ~ l 9 1 /ft , y 7l V /2 . . . ft , y 7 2ny^2n\ 

Wint/C— [2n)\\ 2 / 1 1 Ua 2n&2n 2 a2„ 1 fen / C 

9 



2n \, A 4 y ^a 1 bi,a 2 b2^ a 2b2,o,3b a M a3b3 aibi ---I a2nb2n ai j }1 

q 2 D\ 

TrM n . (B.6) 



2n A 



Here (• • • ) c denotes the connected diagram. Tr'M" has been calculated in 
(1A.17j) . Now we can sum over n and obtain the effective action, 



iVA 4 (A^ 2 — 1) / A 4 \ A^ 2 — 1 / g 2 N 
S/D = --—+ log l— T + — — log 1 + ■ 



8g 2 4 to \g 2 N J 2D * \ A 4 

N 2 (N + 1)(N -3) °- 2 ^ 



Af 2 (N- l)(JV + 3) 

~8~D "°\" A 4 / ' ~ VD ! 



+ „ ^- ; »^^ log(l _^_ +Q (R7) 



C Evaluation of a toy integral using a com- 
plex saddle point 

In Section [31 we evaluated the partition function (13. 8p by a saddle point 
method. A similar calculation was done in Section HJ In this Appendix, we 
illustrate the procedure by considering a toy example. 
Let us consider the integral 

dye~ ay2 - by4 (a,b>0). (C.l) 



25 See footnote [24] for another motivation for denning this propagator. 
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This integral can be evaluated by using the Bessel functions. Alternatively, 
if 6 is small we can expand in powers of 6 and obtain, 



7T 3by/n 
a ' 4a 5 / 2 



(C.2) 



Now we try to solve this integral by using the auxiliary variable. We can 
rewrite the integral as 



1 f * * ( 2 X i o • 2 



dydxexp —ay — + 2ixy' 

Tib 7_oo V 6 

1 f°° ( x 2 1 \ 

= \A 7 ^ \ — b ~ 2 l ° g ^ a ~ 2iX y ' < ^ C ' 3 ' ) 

Let us try to evaluate this integral by using the saddle point method, in the 
limit 6—7-0. The exponent has two saddle points x = —ia/A ± i\J a 2 + 46/4. 
Since in the limit 6 = 0, the extremum is at x = 0, we should choose the 
saddle point corresponding to the "+" sign. We get 



reproducing the earlier expression. 



D The F-Propagator for d = 1 

In this section, we derive the Y 1 propagator in the d — 1 model. The kinetic 
term of Y 1 in (12.151) can be written as, 

r*Tri[y' (-D 2 + A 2 ) Y'] 

= f £ ^ - ^-"l + (a, - «,)• + A') y',„ (D.l) 

where we have used the constant diagonal gauge A 0i j = and we have 
expanded Y 1 ^) = Y^, n Y^e 2mn ^ . We have also used the notation B = iA 2 
as in Section HI Hence the propagator for each mode is given by 

1 1 

(YniAl) = 4.n(a J -a < ) , ]2 A 2 <Wm,0. (D.2) 

— a5 jg r l,a 3 - — aiij + 
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Then the propagator for Y 1 ^) becomes 



<y£w£(o)) = £ 

n 

= E 



e 



2tt / /j j. 



SuS jk 5 IJ 



4vrA 



e 



i(Q!j— Qi)||t| 



Z7T Z7T 

„-A||t|| 



2A 



2tt + Z "2tt 
„A||t|| 



5u5 jk 5 IJ . 
(D.3) 



Here denotes \\t + n/3|| = t for < t < (3. In order to derive it, we have 
used the formulae! 



E 



sin a — nix 



a — n 



71. 



E 



cos a — n]x 



a — n 



7rcot(7ra). (D.4) 



We can write the expression (ID. 31) further as 
M(f)3tf(0)> = 

1 oo oo 

_L e (i(aj- OiJ-AJIItliy^^jyt +e (-i(«j-o i )-A)(^-||t||)y^ a .r» tt ^ 

^^^^ 7i 71 ^^^^ Tl 71 



n=0 



n=0 



(D.5) 



where x = e ^ A and u l n = e tl3nOLi which satisfies Y2iLi u n = Nu n . 



E All loop corrections up to 1/D in the d = 1 
model 

In this appendix, we will show the derivation of the effective action of the 
d — 1 model including the leading 1/D correction ( 14.211) . This correction 
corresponds to iS^o in (13.131) in the dimensional model. Even in d — 1 
model (I2.15p . the same diagrams as in the dimensional model (Figure [8]) 

26 Eqns. flllljl, (|E^5|) and (|R33| arc shown in [39]. 
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will contribute. Therefore, as we discussed in Appendix |Bl the following 
composite propagator is convenient to calculate the loops, 

£ PiW£V))<yf k W£V)) = DNj2G%e^^6 tq 5 kl . (E.l) 

j,p,I,J n 

Note that in this propagator we have only taken into account contractions 
which corresponds to planar diagrams. It turns out that the effective action 
obtained in this way corresponds to the leading term in a 1/N expansion. We 
will make a brief remark about non-planar terms at the end of this Appendix. 

We can calculate this composite propagator by using ( 1D.5I) and a formula 
for Fourier integrals involving 

— / / dtdt'e^-t^e-— l e-— 1 = 5 n+m . , . E.2 

f3 2 J J s/3- 2mn 

Then the composite propagator can be obtained as 

^ijfc = g£2 (^n,ik^ik + Pn,ik^tk + Qn,ikSQ,ik) , (E.3) 

where the n-independent quantities are given by 

oo 

5 « = 1 + Xm ( U -mUm + U k m U. m ), (E.4) 
m=l 

oo 

$tk = (S~ k y = 1 + J2 xm « u -m + u k _ m u m ), (E.5) 

m=l 

oo 

S Q , lk = xJ2 X^iui+m+luUutrniuU - u\) + U- { i +m+1) u\u^(u k 1 - u\)], 
l,m=0 

(E.6) 

and the n-dependent quantities are given by 

1 -1 ^ 1 1 

P~, = — ^-r^ — , P+, 



n,ik ■ i/3(a k -a z )-2Al3 _ ' n,ik ^ i(3(a k - ai )+2A(S _ ' 

27TJ 
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By using the composite propagator (lE.lj) . we can calculate the loop cor- 
rection to the effective action as we studied in Appendix [B] The (n+ l)-loop 
correction in Figure |8] is given by 

( \n 00 N , x n 

-d»U- (fy'DN)" Y, E ( G £«) • A 5 - 1 *) 

m=— 00 i,j=l 

where <i n is a factor derived from the number of the planar diagrams and we 
can fix it by using (1A.15I) . 

di = — 1, d 2 = 3, d n — 1 (n > 3). 

It is difficult to evaluate (1E.8P in general. However, we are interested in the 
theory around the critical temperature x ~ where x is given by (14. 9|) . 
Hence we can expand the effective potential with respect to x and the lowest 
order of x is enough to evaluate the l/D correction. Especially, only the 
coefficient of x|«i| 2 and x 2 |mi| 4 and the gauge-field independent terms will 
give us the relevant information of the dynamics around the critical points. 

E.l Two-loop correction 

The two-loop correction to the effective action corresponds to the n = 1 term 
in (1E.8|) and is given by 

N 00 

-\^dny: e g %- ( e - 9 ) 

i,j=l m=—oo 

We sum over the Fourier mode first. 

00 

1 + 2 J2 ^ m <y- m , 

m=l 

(E.10) 



E P ^ = 1 + 2 E 



a -m a mi 



m=— 00 
00 

E Q 



ni = l 



u l _ x + U 3 _i 
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where we have used f]D.4|) . After summing over we obtain the correction 
to the two-loop effective action as, 

^ + 2^(x 2m + 2x m )| Mm | 2 j +S int , (E.ll) 

(x 2 (u 2 1 u^ + u 2 _ 1 u 2 ) +0(x 3 )) . (E.12) 

Here we have omitted higher x terms in the interaction, which are irrelevant 
around the critical points, as argued before. 

We notice that this correction includes a cubic interaction x 2 u\u_2- Since 
the effective action (14. 6 \ has the term |u 2 | 2 /2 arising from the gauge fixing, 
| Mi | 4 term can be induced through this interaction after integrating out u 2 . 
However, the coefficient of the |wi| 4 terms obtained this way will be 0(x 4 ) 
and we can ignore it compared to the \ui\ 4 potential in (I4.26p . Generally, we 
can show that the lowest order coefficient of the cubic interaction from the 
higher loops is also x 2 and we will ignore them here. 



gtwo—loop 



N 2 /3A ( A 



S 



N 2 (3A f A 



hit 



4 \ A 3 



E.2 Three-loop correction 

We evaluate the three-loop correction (the n = 2 term in ( IE. 81) ) 

m=— oo i,j=l 

In order to calculate the sum of the Fourier mode m, we derive the product 
of P ± and Q as 

i 

P~ . .p + . . = (p~ . . + p+ ..) 

m,ij m,ij 2A/3 m '' l 3 m i*J/ ' 

Pm,ijQm,ij = {P m ,ij + Qm,ij) , P m ,ijQm,ij = ^~P m ,ij + Qm,ij) ■ 

(E.14) 

We also calculate the sum of squares of these quantities by using the formula 

°°/i\2 cxd , 2 

E (—) =-1- E — = < E15 > 

\a — m J oa ' a — m sin na 

m=—oo x ' m=— oo 
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This leads to 



n= — oc 
oo 



'1 — u^^x 2 ) 2 



4x 2 M i 1 ni 1 



4w!_ 1 iq 



(1 — n'^ui 1 x 2 ) 2 

(E.16) 



(1 -u* X) 2 

m=— oo v — -L 1/ 

Now we can sum over i, j and obtain the leading order of the corrections as 



gthree—loop j\j2 



3 (5\ 2 
128 A 5 



+ 



32 A^ 



x + 0(x 2 ) |tti| 



3 /3A 2 2 /5 



32 A 5 



+ /3A +0(x 3 ) | Ul 



(E.17) 



E.3 (n + l)-loop correction to effective potential 

Up to three loops, the leading order of the coefficient of |wi| 2 is x and |ui| 4 is 
x 2 in the x expansion. We can find that these are true even in an arbitrary 
loop. Thus it is enough to fix the coefficient of these terms in each loop. In 
order to evaluate the (n + l)-loop, we have to calculate 

oo JV n , ^ x 2 oo N 

E E (pm,ij) = (g/^J E ^ ( P m,ii S ii + P m4i S ii + Qm,ijSQ,ij) 
m=— 00 i,j=l ' m=—oo i,j=l 

(E.18) 

in (IE.8p . However we can reduce this equation. Since Sq^ is an order x 
quantity as in (IE.6j) . we can ignore Sq (k > 3) terms here. Then we should 
evaluate only 

(ps- + ps + ) n + nQs Q (ps- + ps+y 1 



+^-^-(Qs Q f(ps-+ps+y 



(E.19) 



where we have omitted the indices. 

First we calculate the first term in ( IE. 191) . It is convenient to define 



o-k,i = (P + ) k (P ) l so that the equation becomes, 



(PS' + PS + ) n = (f)a n - k , k Sl- k S k _. 

k=o ^ ' 



(E.20) 
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Through the relation (IE. 14[) . a fci ; satisfy 

a k ,i = (ak,i-i + Ofc-i,/) . (E.21) 

Here we can approximate a^o = ao,fc = if k > 2. This is because we 
can show that they are order x 2 quantities and the lowest order terms only 
contribute to x 2 |mi| 2 by using a similar logic as in (1E.15j) and (IE. 16[) . Then, 
through (1E.21I) . we can obtain 



k fc3 &2 



V ' ' fc,_i = lfc ; _ 2 = l fcj=lfci=l 

(ai l0 + M • (E.22) 



(2A/3) fc +'~ 1 (fc — 1)!(Z — 1)! 

(lET20|) as 



Since ai j0 = -Pmij anc ^ °o,i = P m iji we can sum over ^ ne anc ^ m indices in 



oo N 

. n—k 



m=—oo i,j=l 

=2N 2 (1 + 2nx\u 1 \ 2 + n{n - 1)x 2 |mi| 4 + • " ) » ( K23 ) 
where • • • denotes the irrelevant higher order terms. Then ( 1E.20|) becomes 

(1^1, g I'll; (1 + 2H".I 2 + " D*>.l 4 ) + ' ' ' • (E.24) 



where we have used 



fc=0 

Next we evaluate 

,n-l 



QS Q (PS- + P + S + ) n \ (E.26) 
in (IE. 19[) . First we can show that 



oo N 

/] / ] (Pm,ij) {^ij) i^ij) &Q,iji 



m=—oo i,j=l 



42 



does not contribute to the relevant potential. Hence, by using (1E.14[) . we 
obtain, 

oo N 

^ ^ Qm,ijSQ t ij (P mt ijS i:j + P^ijSfj) 
m=—oo i,j=l 

j oo jV 



(A/3)' 

v ' m=— oo i,j=l 

4iV 2 2' 



2on-l 



Finally we evaluate 

(Q5q) 2 (P-S'- + P + S+) n ~ 2 . (E.28) 



Here we can show that 

oo N 

£ E(4)K)'(s;)'^ 



m=— oo i,j=l 

does not contribute to the relevant potential. Thus we obtain 

oo V 

^ ] ^ ] {Qm,ij£>Q,ij) {-P m ,ijSij + Pm,ijSij) 
m=—oo i,j=l 

n _ 2 oo AT 2 

\ ' m=—ooi,j=l ^ ' 

(E.29) 

Let us summarize all relevant terms of the (n + l)-loop effective action. 
The gauge-field independent term becomes 

Note that this result is exact. The leading \u\\ 2 potential in the x expansion 
is given by 

-(-Td n N^A (^-) (^m + i y M + 0^). (E.31) 
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The leading \ui\ 4 potential is 

-(-) WA^l (A)" (I^H + 2 + A,) + 0<*«). 

(E.32) 

We can sum over n by using the following formula: 

^H^^^^v^-l, (E.33) 



71=1 



and its derivative with respect to x. With this, we finally obtain the effective 
action (1QT]1 . 

Note that it is possible to extend the calculation in this section to finite iV 
case as we have done in the d = model. The finite N result for the gauge- 
field independent constant term is simply obtained by replacing d n N n+2 in 
(IE.30j) with TrM n in (IA.17j) . The terms including the gauge potential are 
more complicated and we have to modify the composite propagator (lE.lj) . 
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